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fA  /i  I 


utils 


obit 


The  motivation  for  this  study  is  a recent  pa:  :r  of  "iff,  Fchcr.r.,  e 

and  Varga  [7]  where  it  is  shown  that  there  exists  a sequence  , 

with  R (x)  - P (x)/(l  + —)?  P , Z H such  that 
m m-1  rr.  m-1  m-1 

3 - 2/2  <~;'e-X  - R |!  1/n  < i 

^ L [0,-) 

That  is,  {R  (x)}  ccnver  . geometricallv  to  e " on  £3,  ®)«  In  addition, 

since  the  soles  of  R ( ■: ) a;  e all  real  it  follows  that  R :z)  must  c: eve:  : 

n m 

geometrically  to  e “ in  an  infinite  sector  symmetric  about  the  Positive 
x-axis  [S].  An  - li cation  of  this  theory  is  in  the  construction  cf 


numerical  solutions  for  solving  linear  systems  •'  rcinary  Jif ferential 
equat ion 3 whi  : jr ise  free:  : ni-discretization  r f lin  ar  paral  : z partial 


differential  neat  ions  (see  Cl]  and  C 7 3 )-  For  ex  imp ie,  ar  iescr i 


: t-  J r-  r 7 ’ 


consider  the  n • srical  solution  of  the  linear  syst  m f ordinary  ii  : 
tial  equation.' 


(1.3) 


du(  t ) 


- -Au(t ) + k,  t > 0 


u(0)  = 

-0 


where  u(t)  = [u,(t),  ...,  u^ft)]1  is  a column 


ve  o tor  with  i c it  nent 


and  A is  an  hxn)  positive  finite  symmetric  matrix.  The  integer  n is 
related  to  the  m<  si  air,  of  1 ii  cretization  and  can  be  larg  . The 
solution  to  (1.3)  is  ~ ever.  • ...  licit’;/  ” 


(1.4) 


i(t)  = A.  ‘ k + exp(-tA)(: 


for  all  t 0,  where  cx.(-tA)  = l (-t.A  ■ /v!.  Fcr  cc~;  utaticnal 

v=c 

one  m jst  appr  oximati  •••  In  [7]  I • by  using 


fL(dtA)  j (I  t — A ) np  (AtA)  wh  re  i € It  , is 
m f.  m-1  m-i  n-1 


>•  ?o;ut 


. cf  II  -x  P(x) 
inf:  he  .. — . 


(1  + -p)“  L C 0 ,~ ) 
m 


: r e n ,}  = 

m-1 


TIi*  - •;  sue;,  a ? 


4 


Which  Will  involve  increased  computation , where  R*(x)  = ?*  <>;)/ 


m 


x ) . 


We 

say  ap 

R* 

e . 

m 

in 

as 

corres; 

-3  * T — • - 


i=I 

ir.ce  our  nur.or ic3. 1 results  g ~ m ~ 

ria-  13 * r-£  appears  to  give  a rise  to  the  come  sort  or  r.ec'r  *d 

reasec.  accuracy  for  no  additional  ef for  - . • 

this  is  known  to  he  true  in  theory  also,  for  t - case  that  m=2  [4]  (R.  S.  Varga  has  inform*  : ; 

this  has  a.  s . ; ■ r 'r  ■ - : Wi  * ..  . „ \ ..  .. , , 

" ' * '•  we  will  return  to  this  ise  iter. 

In  the  next  tv;;  sections  s'al1 
characterisation  ar.d  unicueness  ^.,-.^>3 
the  special  o ->3 p 


tenc‘ 


note  and  > . -•  * -cr  - 

•n  - '1  ’ * c’~ 


.se  c.  atprc.xi~.anrs  •->  i'"1  <*=1  rr--«  ^ -v--  ...  ..  -- 

->  - •*  * U » / -I-..  . . . - . ;<2n  V/6 


describe 


Sll’  r -2^-4 


numerical  ro_i 


some  epen  pro'  lore. 


- tn«  j.ly,  we  .;i  1 close  th  ; : it  »r  with  a 


2.  Theoreti  / : 


Vie  be. "in  t 


from  f 


1 tms  section 


■‘1  th  a 7 reel  cf  existence  of  b •; 


1-  i r h f e 


f £ C[0,  a]. 


Th co n • ~.  . . . . ; ix  1 C ^ ^ — 1 ~ — v -a  **  -p  • .„ 


**!!  --  i-f{  ;|f  - : K % } wh-  re  , 


1 . li 
11 


II . 

•1  11 


Proof:  Let  us  a3-.ee  that  f <£  f y 0 c V 
(i.e.,  j|f  ;f  > Wf))  and  t .it  n s:<{:  (x)  : 
this  last  cor.dit  i - i 


is  not  a best 
x €.  [0,  01 ] } = 


L [0,ct] 

approximation  . : 

•j  r J » %•  .+. 


' ” -S 


— rj  ^ r>  *- 


below . T'e; s , 

[a,  b]  C (C, 

T 6 + b 

e"  xo  = T* 
Y ♦ x 


‘1 J 


1 - — 


- s:<-  ;“3  3 closed  interval  [a,  b],  b ' 

: £.  [a,  h]> 


i ~ in{ f(x) 


a,  suer,  tr.at 

= Y > \_  ( f > . 

th.it 


IS 


(A  = \ ( f ) ) . 


(2.1) 


Y + ^ y + X. 

:(x)  - — —— 1 + — o 


for  all  x € [0,  a].  Now,  let  ua  normalize  R.  (x)  by  requiring 


Q,  (x)  = H (qW(x  - x„)  + 1)  where  0 
k i=l  1 ° 


< 

— *1 


(k) 


Note  this  can 


done  since  Q^(x)  is  known  no  have  all  negative  roots.  Thus,  if 

rj 

Q (x)  = II  (x  - r.),  o <_  m and  r.  <0  for  all  i = 1,  ...,  p then  we  may 
k i=l  1 1 

rewrite  it  as 


Q.(x)  = ( n ( x_  - r.))  n [(— — )(X  - X ) + 1]. 


i=l 


i=l  *0  ri 


Sot  q<!'>  = — - — 
i x ,-r. 

0 i 


to: 


(k) 


« t-1  > Q ^ 


0 fcr  i = p + i 


(k) 


and  note  that  r.  <0  implies  that  0 < c'."'  < 

l ■ 


: or  i 


to  he  incorporated  into  F.  (x). 

K 


Finally,  the  constant  ( !■  - r.)) 

i = l v 

Since  f is  bounded,  v:  have  from  (2.1)  that  there  exists  a constant 
M > 0 independent  of  k such  th-t  for  all  x £ [0,  a] 


(2.2) 


-M  <_  F^(x) 


. . Ik),  m 

Since  l c . ; . . . _ , 

* 1 1 - -L  , . 


C [0,  — ),  we  may  xtraot  convergent  subsequent - 
X0 

(k) 


(relabelling) such  that  q.'  ' -►  q.  to  [2,  — ] fcr  i = 1,  . . . , m.  Note  f 


li  x. 


then  a.(x  - xn)  + 
-i  u 


•*  -i  - r 


to 


Thus,  Q.(x)  convert 


to  Q on  compact  subsets  of  the  real  line.  Now,  (2.1)  restricted  to  [a,  1] 
gives  that  there  exist  constants  both  positive  and  independent  cf 

k such  that 


(2.3) 


c.Q,  (x)  < F..(x)  < o - Q,  ( x ) 


. Cr 


for  ail  x £ [a,b].  New,  for  x £ Ca,b]  at. ' 


CO,—]  vo  have  th 
:'o 


above  proof.  This  is  so  since  in  the  case  lit  f(x'  - i'c'i  and 
| f ( x ) j < ||f|j  for  all  x G [0,  then  for  r,  sufficiently  larpe 

the  interval  [r.,  n + 1]  car.  be  used  for  the  interval  [a,  b]  provided 
0 is  not  a best  approximation  to  f. 

Now,  we  wish  to  study  the  space  , In  what  follows  we  shall 
prove  a local  characterization  and  local  uniqueness  tnocreo  for  this 
space. 


Definition  2.3.  ?.(:<)  r (p^  + F2X  + • • • + Pm>:  )/ (qy-  + D /l-ra  c_i  i; 


best  approxirv. tier. 

that 


€ C[0,  a]  £n  [0,  a]  is 

m-1 ,,  - 


R(x)  = t p x * ...  + prx  >*.  qx  ■*  ire  anjj.  \ , - ? L\  < 6, 


i = 1,  . . • , m and  |q 


£ phen  ||f(x)-R(x)  !j  < ||f(x)- ?.(y.) | . In 


addition,  if  strict  ' n - - p . c 1 1 v_  h o : '.a  'whenever  R(x)  f R(x ) then  F 


to  be  lor a 11 v uni  . 


Before  we  can  prove  cur  local  characterization  tr.ccrem,  v.  .must  err. 
two  lemmas.  The  first  lemma  states  that  P has  a local  Kertnl- 
property  of  order  cne. 


Lemma  ? . u ■ Sucre  ce  ?(:■:)  - ?(x)/Q(x)  - (pn  + p x + ...  + p^x  ). (qx 
is  nondecenerate  (i.e.  ’ i 0 and  P and  3 nav  nc  com:  r.  fac  s) 


Let  0 < m,  <_  ra  and  m.  + r„  = m + 1.  Sutpcr  ' {y.  C [0,  aj  with 

yj  < yXll  f-r  ~2--  ihl 

5 > 0 and  6 > C ruth  • hat  if  1 y . - y . j <_  5,  and  In.  - R(Cr)|  < r.  for 
i = 1,  ....  m„,  [y;  . - yu  ! < 5,  am!  \z\  - R'(h,)\  < ' fyr  j = 1,  . 

^ j J ”*■  J j 

then  there  xi  '-t-t  L" 


R(x)  * ?(xVQ(x)  = (p.+rvx*...-„xm""Vf 

-i-  *- 


with  |p  - p i < 5,  v * i,  ...,  m,  ;q  - c|  < i,  R(y. ) = z ; , 1=1.. 

— J V — “ i 


■ 


R’(y[.)  = z_:,  -• 


p.  , . • . , , q dc  en tit 


8 


Proof : This  result  follows  from  an  application  of  the  Implicit  Function 


Theorem.  Thus,  one  forms  t 

3 - (p, , • • • , P » q,  '/,  > • . 


.e  system  f fa)  = 0,  u = 1 

y " 


m + 1 whei 


“ T I XT1 

- t.o  » / 1 » •••  j y "ht  $ 


*-■’»  • • • » /-r 


tr  * tr  1 

TT;o  » i » • • • * *-•; 


’ 1 


•T-  1 


-U'-'  a *Vu  " (qyy  + irzu’  p = 1»  •••»  ^2  and 


f„(a)  = p,  + . . . + o y’“ 
y - 1 ‘ nr  v 

^ir.o+y^J^  * ' + D(?t  + ...  + (ir.-l)p  (y')  ) - rr,q(p1  + , . . o (y ' ) ") 

jl  u nip  i ‘nyp 

y = 1,  ...,  m, . Observe  that  the  point 

” • r ( v ) RCv  i v ' ( ~ . 

" -imp’  -r  ’ •••’  "'J^2  ’’ 


- (qv’  + l)m+1z’, 
y y 

aA  = ‘ 


.9  ‘l  9 J •,  5 


' m 


2 "1 


R'f’/i  ) ) sat  it  ties  this  svst . 

ml 


- hus , since  eacn  r unction 


: n i s 


system  has  continuous  first  partials  with  rc-spe 
of  a (or  components ) v: . need  enlv  move  th.it  J a 
with  respect  to  p , ...,  p , a has 


cn  or 

t he  va 

J(a), 

cf  the 

a nonzero  cv "remnant  at  a - . 


by  using  the  equalities  z.  = ?,(y.- ) and  c!  = R'(y.-  ) and  a^diro-  (-.  * ] )" 

3 j 


times  the  v'h  row  to  the  yth  row  where  y > 
becomes 


and  y 


Vi-"  2 


_ra.-l 

yl 


v 

- m2 


det  (J(a  )) 


?(y)  > 


q qy:  .-Key;  +1) 


q qyi_  Kqy;  *i) 

mi  *ni 


— m-l  . , 

(*• 


+ _ 1 


-mv  r-- 

■'Kqy.  +x ) 

i 

P(yn, » 

■"ytn-.TT?  ZTT 

* ' u rr.2  - 

-V'-V- 


iy:  +(m-l)y-  ' ( .y-  +1)  -n~:  ?’(“•  )- 

"n]_  -tip  *rip  ' iI7li  ' lriq 


Assuming  that 


to  the  resulting  function 
Define  :!(":)  by  H.(t) 


•1  > C,  replace  in  the  *.  +Ist  r 
‘*‘1  / 

Mote  that  C C-  v. 


det( J(e  ) ) 


"«•  by  t and  set  G(t)  ecu 
..  and  ) is 

row  ret. laced  b ’ 


L 


9 


( ( qt  + 1),  t(qt  + 1),  t‘  (qt  + 1),  -mtP(t)).  Note  that  H s.  ii^, 

H’(t)  = G(t)  and  H(y  ) = C , v = 1,  . . . , m, , «'(>-.)  = 0,  j = 2,  ....  n 

V 4l  j -i- 

so  that  H has  t.  zeros  counting  multiplicities.  Thus,  if  H t 0 then  H 
can  have  no  more  zeros.  Hence,  if  we  can  show  H t 0 then  it  will 
follow  that  H'(y?,  ) r 0 and  so  det(J(a.-,) ) i 9 as  desired.  Now 
H(-  n.)  = (-1 )‘"+  "*■  -_-F(-  -1_)3  where  D is  the  determinant  obtained  from 

q q q 

s +-  st  . - 

det(J(a„))  by  deleting  the  m+i  " column  ar.d  the  m.,tl  rev.  Since  r an. 

- 1 . 

Q have  no  common  factors  we  have  that  ?(-  --)  i C.  But  now  aiding  (-c) 

q 

tb  th  _ 

times  the  v row  to  the  p~‘  row  where  p > m and  v = v:  (note  the 

2.  j ■‘•p -mo 


row  contair.it ■' 


y;,  is  gone  rrom  P , this 


yi_,  ...,  y ) and  then  factoring  cut  ( qy  ,•  +1)  free,  rtw  m-  -r  p - 1, 

p = ?,  . . . , nu  shows  that  D equals  a nonzero  constant  times  a detei  ni: 
which  is  known  to  have  a nonzero  value.  A similar  proof  v orks  for  the 
case  that  = 0 (no  derivativ<  - present).  In  this  case  cue  simply 
replaces  y ir.  the  first  row  by  t and  proceeds  as  ab  :v«.  without  refers*  I 


to  derivatives.  Finally,  to  guarantee  that  R £ y we 


outre  tr.at 


5 < |q|. 

Lemma  2.4  gives  a tointwise  local  solvency  preyrty  when  m,  = 0, 
pointwise  in  the  sense  that  the  and  £ depend  upon  the  points  at  whi 

the  functions  are  being  evaluated.  In  order  to  prove  the  necessity  of 
our  local  characterization  we  need  the  following  zor' -count ! r.t  trotert; 
Here  we  shall  assume  that  a > 0 is  finite,  a < a. 


Lem  a 2.5.  Let  ?<:<)  = ?(x)/Q(x)  = (px  t . . . + pmxn,-A  )/vqx  + iff.  £. 
be  nonde  gene  rat  and  q > 3.  Sucrose  0 < y,  < < ...  < y < a and 


y d to,  n] 


>01. 


10 


to  P/Q  and  the  y-'int  set  (y,  , . . . , y , y } accord:--  to  Lemma  C.u  vi -h 

m = 0.  F Ira] l.y , for  all  a,  jaj  < 5-,  let  R (x)  = ? (x)/Q  (x> 
x x o a a 

= (pla  + •••  + Pma5<m  + ir  — lii£.  e ^ liJnJi 

ga^sfi^  R (y  ) r R(y  ),  i = i,  ....  m,  R (y)  = R(y)  + a,  jp.  - p.  J < {. 

u i a ^ c i — 

i = 1 , . • • , rn  and  ' q_  - q j <_  5 . Then  there  exists  <$„  > 0 such  that  if 
0 < |a|  < 6.,  then  the  cr.lv  z—  of  ?.„  - R in  [0,  a]  are  y,  , . . . , y 
and  — R changes  sicn  a v.  each  of  th/sc  that  sra  in  (0,  ci), 

Rroo  i ' oiip  l OS6  not.  > p.on  t hor  e exists  ^ 0 o • 0 f on  i suc1?i 

3 3 


that  R, 

0 . 

- R oith 

er  has  an  additional  zero  at  vg  G [0, 

x ty2 , 

• • • 1 y_ j 

3 

3 

-a 

or  R3_  - 
3 

R fails 

to  change  sign  a4-  one  of  the  points  y 

V 

£ (0,  a). 

In  both 

cases  we 

write  yg 

for  the  additional  zero  with  the  und 

3 

t . . Ling 

t r.at 

Vo  . = y ; 

3 * 

for  some 

1 means  that  ? _ - R dees  not  change 

sign  at  y 

- "n  ’ c- 

case . By 

passing 

to  a subsequence , we  may  assume  that 

Yg.  + y * £ 
j 

[0,  a] 

where  her 

e v;e  are 

using  cur  assumption  4 hat  a is  fir. its 

. We  new  c 

oncider 

two  cases . 

CASE  1.  y*  £ {y.;.:sl.  Choose  5*  > 0 and  > C corresponding  to  ?/Q 
and  the  point  set  / y - , . . . ^ ’/'*}  seconding  to  Lo~'  a i w**’1,  rr' , « q 

Then  for  j sufficiently  large  we  have  that  jy.  - y*j  < £;•: ^ 

l5°]<i'0j>  ' I;i';  • ^ li5'  f0r  1 * 1 ” ^ 

|qa.  “ <.5*  tine:  the  parameters  p.  , ....  p , c depend  continuously 

' j.<3  rr.  a 

upon  the  refraining  parameters  (so  **:.  it  ?.  converges  uniformly  ? on 

* 1 

compact  subsets  o:  [0,  ■; ) ) . Mow,  Rj.  and  R agree  at  the  points  y,  . 

so  that  by  the  uniqueness  part  of  Lemma  2.4  we  must  have  P-  = R which  i* 

j 

a contradiction  since  F,  (y)  t R(y). 

Tj 

CASE  2.  Suppose  y*  - y,  for  s me  J. . Choose  <5{:  > 0 an  - 0 > 0 -nr"  ---• 


tc  the  point  set  (v y } iinc  to  I 

* • ^ ’ 7 p 


1 


e that  R.j4  - P. f 

j 


Then  for  j sufficiently  large  we  have  that  jy0.  - y.j  <_  c*. 

J K ■) 

for  all  y in  the  closed  interval  It  with  endpoints  and  y . , u - 

i = 1,  m and  |q0.  - q|  £ o'--.  Now  from  the  fact  that  - R’  is 

continuous  on  I.  and  P0 . - ?.  vanishes  at  yc . and  y we  Lav-.  by  F.oiuo * 
j ] J *• 

Theorem  that  vanishes  at  some  point  % • . £ Ij  provided  yai  * y 

If  va.  = y.  (for  some  j)  then  y £ (0,  a)  and  we  have  that  f ,L  - is 

1 *.  *•  J 

zero  at  y„  since  R~.  - R does  not  change  sign  at  this  point  in  this  u 
Thus,  yj  £ I.  with  P-Oj  j ) = -n  either  case.  Also,  jy’,  - v 

and  |R£.<yr.)  - R'(y,)|  = | F.  ’ (yj  . ) - ?’(y  )|  <6*  Thus,  by  the  uni 
part  of  Lemma  2.-  we  must  have  that  R0_.  = 7.  for  j sufficiently  large  s 

J 

these  functions  agree  al  v.  , •••,  y and  their  derivative:  : e at  y„' 

which  is  our  desired  final  contradiction. 

With  these  result  ■ w are  now  ready  to  prove  our  local  characte  i 
theorem  which  is  an  alternat ion- type  r-*ault . 


- V < 
2.1- 


Theorem  2.5.  L ■ t n > C.  I 
with  q > 0 is  a local  ”_yot_ 
if  an  d or.lv  j y ^ __  ■ rv.r  c_. 
altern  it : " t extr  me  noir.ts 


, _ _ m-1 

nen  a r i • ■ . • 1 ••  ’ x v 

approximation  to  f € CCO,  i]  on  [0,  a]  fro: 

:rve  E(x)  = fix)  - ’’lx)  has  t !c y n + I 

. ( If  a = ®,  then  w ■ r cuire  ' i-"  f (x)  - 


Proof:  The  nec  .ssity  of  this  alternation  now  follows  by  tr. : ary-:’  n * o 


of  Thcore::  7.3  [h 

, pp.  lo- 

12] for 

v ar i s lv  r.t 

' ''  nS  <r  ' 1 : 

number  correspond! 

ng  to  the 

decree 

L- 

re),  2 • 5 • 

above,  is  need-'d  f 

.ir  cons  tr 

uct ing 

a bettor  a 

has  less  than  m + 

2 aiterna 

•:  in  ; x 

Lnt 

. If  0 and  a are  both 

e x o : ;me  jint  . 

v t r ; f 

ervuri 

: n ensi  i a 

.5  may  1 n<  : 

for  the  case  tnat 

a - ri  ^ w'-’ 

r.:t  i t 

hat  since  . . ' 

Vi  ( f ) 

went  y rep]  [0, 

h-  [0,a.l, 

i ' ’ 

- for 

x>‘,  ( - 
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which  is  our  desired  contradiction.  Indeed,  there-  are  two  ibi.ltio. 

to  be  considered. 


CASE  1.  y*  < y*  < ...  < y:’:  . In  this  case  we  apply  L rein  2.-  wi*h 

= 0 to  P/Q  with  respect  to  these  points  (i.e.,  y . of  Lemma  2.4  i y£) 

with  y.  of  Lemma  2.4  set  equal  to  y..,,  i = 1,  . ..,  m t 1,  ' fix  »:.d 

4 1 J 

z.  = R(y..),  i = 1,  ....  m + 1.  Then,  for  j sufficiently  1 .rv-  w. 

1 j 

that  R.  satisfie  the  : nclusioi  2 . M (i.  . , eff  i i f • 

3 

sufficiently  close  to  rest  ?ctiv  coefficients  of  R r.:d  = / . with 


iyii 


•vl  and  |z.  - R(y$)  small).  But  R also  sat isf i t : >n  3 . : - 

1 1 ' i - i ■ 


of  Lemma  2.4  a:.!  since  botl  ' ' wi  hav  by  • - uniqu  n :i 

m 

Lemma  2.*  that  R.  ‘ R. 

J 


CASE  2. 


car  v: 


-'ni  " ■'n1+l’  • • ‘ ’ -'nj.  - “ 

with  ns  - £ to  the  points  {y.  , ...,  v } (a  litti:  • . tl._  11- tin  t 

i -1.  m-  it  i. 


points 


1 


• » an  : th  j int:  {r  .,^ , 


Le  1 1 1 r.  * u . 1 


the  first  ind-.-x  v sue!  that  y -*■  y.  as  k ■-  «,  i = 1,  m - l + 1, 

we  take  the  y.  of  Lemma  2.4  as  y , (k  fixed)  = R(y  . ), 

-1  .•  t.  i • p s K 


i = 1 , . . . , m - £ t i . We  also  t?V- 


'••S* 


definition  "ust  mot  to  c:i  > . of  ;ur  : roc:  .r  . r. . = R (v  • 1, 

‘ i 111  K 

j 


j = 1,  ...,  i.  The  desin  i re  lit  th  n : 
with  k playing  the  role  tl  j lid  in  la  • ! . 


a r.  ; r.  . 


Corollary  2.7.  Stir-.  ? > l’.(x)  - Ti  : (:•:)  = ( 


h — 1 


..  t prx  i‘)/(qx  + l)rv 


is  a local  lest  a. 


f(x)  from  ( n [0,  o]  • ■ 


wi th  q > 0.  Then,  R is  I nr: 


The  converse  of  Corollary  2.S  Is  most  likely  false  since  in  the 

“X 

m = 3 case,  f(x>  = e and  ct  = ® we  have  essentially  at  least  two  b t 
local  approximations:  R^(x)  - (1.0CSC5  - .2701Cx  + .01447x‘ )/( .27127:<  + 1) 

-3 

with  error  norm  3.05002  x 10  , achieved  at  the  extreme  points  0,  .462, 
2.178,  6.876  and  37.260  (with  e*  - R^(x)  < 0 at  0)  and  ?.0(x)  = (.28663 


+ 2.52827X  - .44972x2)/(1.051CSx  + l)3  with 


;rror  norm  1.33720  x 10  ’, 


achieved  at  the  extreme  points  0,  .172,  .3'72,  2,950,  13.225  (with 
eX  - R2(x)  > 0 at  0).  These  approximations  were  computed  over  a 20, 001- 
point  equally  spaced  grid  imposed  on  [0,  40],  Although  coefficients 
given  above  were  rounded,  using  the  actual  coefficients  computed  the 
absolute  errors  at  the  extreme  points  in  each  case  agreed  to  at  least 


15  significant  figure-.  It  seems  very  likely  that  a theoretical  argument 
can  be  given  starting  with  these  two  functions  to  show  that  a-  least 
two  distinct  local  best  ap  r '::i::\;ticr.  . :::i.  t for  this  problem. 

We  can  extend  seme  of  our  results  for  K to  the  other  possible 
configurations  of  the  denominator  of  members  cf  if  . If  rr.,  , n,  > 0 

TTi  J.  it 

and  m1  + ...  + m,  = m,  we  define  , ...,  m,  = {?.  = ?/Q  : ? £ ilm_,  , 
nii  m- 

Q(x)  = (q  x + 1)  ‘...(q  x +■  1)  ' , 0 < q.  < . . . < q } . Although  in 

JL  K,  — i. 

general  we  cannot  expec t exist  nee  of  best  appi  xitnations  fr  m 
^ , ...,  since  the  set  cf  allowable  coefficients  is  not  closed,  we 


Theorem  2.9.  let  l > 0,  r . , ...,  m > 0,  and  n,  + 


£ nonde.ger.erat-  f '(y.)  - (p.  * . 
with  q > 0 Is  a loca: 


l r ar  . c n • 


m j * * * • » Ai 
mt ; +1  alt-  • 


fur-.  1 -r-  “ , in  * -.1 


- p xm”A)/(q.x  - : ) .. .(q.x  t 1)  % £ 

X K 

ion  to  f £ Ct  , ot3  on  [0,  q]  fr 
ve  t x)  • > - R(x)  : 


) i i a-  4 


The  proof  ol  this  theorem  requires  nly  prov 

2.4  for  ^ . This  proof  is  ncre  involved  t! 

ml»  •"  »n2. 

2.4,  but  follows  the  sate  lines;  the  variable  : ow 
((q1t+l)...(q4trl),...>t!n"1(qLt*l)..  .(qtr  + l),  -n. 
...,  -m  t(q.t+l). . . (3,  ,t+l)?(t)). 

As  two  consequences  of  this  result  we  note 
erate  bent  approximation  I.  to  f from  )C  with  al' 
positive  is  such  that  f - R har  only  rn  + 2 alterr. 

then  R £ X!  ; second,  : f i nond  • nei  ate  b<  st  icdi 

m 

'P  has  all  its  ienon  inator  : ff  1 • : positive 

m 

actually  the  unique  best  a;  yrexin  • ■ ;;  : ;rr  ; 

P <c  a , Q € It  , Q > C on  [0,  a]}, 
ns- 1 m 

So  far  we  have  always  cor.strt  i:  . th-  .v. . : : 

in  II  . , but  if  we  replace  R , ' and  : : : , 

m- 1 n- 1 n 

expressions  of  numbers  of  alteri  iti: 


’no  no: 


=ce  d ir. 
.his  a: 


3")'  a o 


he: 


• . i ... 4 

In.rv 


1 1 


.I}  I 


r 
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precisely  the  linearization  of  g(q,  , ....  q_,  x)  via  Tayl.r's  thecr- 
applied  to  the  first  m independent  variables.  The/1  al/'.rith.r.  use. 
the  same  approach;  the  linearised  denominator  for  thi  il  orithm  L; 
qmx ( qx  + l)m  + [(1  - m)qx  + l](qx  * l)n‘  1 = qy  (q,  x)  + : .(",  x). 

If  a is  a large  finite  number  or  a = 00  (in  which  we  c i • we  tens i : 
[0,  a]  instead  of  [0,  »)  for  some  1 irge  finite  a)  , and  we  wish  to  i 
fairly  fine  mesh  in  order  tc  get  m i irate  appr >xin  itioi  ter  CO,  sj, 
then  car  .'Hr ) will  be  large.  Since  this  Leads  tc  a large  and  iif  icult 
linear  programming  pi  bl< m an  i nan  si  ra • ■ problems  in  th<  differ- 

ential correction  algorithm  we  us  id  the  Re:  s-Difc  ilgorit  mi  [ J for 
calculating  the  linearized  minimum.  This  alf  rithm  i]  lies  th< 
differential  correction  algorithm  t cert  in  - ill  of  T chosen 

in  such  a manner  (depending  on  iltematicn)  that  icnverg  nee  t > th< 


solution  on  . occur.’. 


• ' . . . . w _ . . — -J-  r. 


work  sir. or  a standard  alternation  the  ry  has  not  develo;  3 foi 


problem;  however,  in  or 


: ' : . inner  ar.d  outer  ilg  rithms  j v<  rred 


and  we  obtained  approximated  satisfying  theorem  2.9.  Alth  igh  a 
precise  stuly  of  these  algor:  thr.s  remains  to  be  done,  we  cc:  jecture  th 
the  A’  algorithm  will  renverg'  (as  umine  the  converg  nc<  of  the  inner 


the  A algo 

rithm  will  - : 

iterations ) 

if  the  ir.it  i, 

sufficient! 

y good.  We  rr„ 

as  present! 

y : t i t ul 

4»i(q1 

q_,»  x)  = o. 

next  sta ; 

will  : * be 

algorithm  1. 

•;  :ho s*  r q. 

*i 

(or  = 0 

if  i s i ) or 

^ i.  ,.  * c 


. then 


algorithm  1 : chose-r  q.  and  , at  the  next  stipe  jo  th.it  eirh  m q.  = 
(or  qf  = 0 if  i = a ) or  - . = i.  (cr  a.  . = 5 if  i = r.  - i). 

*i  -i*. 

As  an  ex  ' . i • . rob 1 \ > t i n . 

p + ' ,x 

: . . ’ ft!  r— — - - - - • wit 


f 


VUI 


